hep-ph/9903439) 
DTP/99/28 
March 1999 



Small x QCD effects in DIS with a forward jet or a forward 7r°. 



J. Kwieciriski a ' b , A.D. Martin" and J.J. Outhwaite a . 



a Department of Physics, University of Durham, Durham, DH1 3LE, UK. 
H. Niewodniczanski Institute of Nuclear Physics, ul. Radzikowskiego 152, 31-342 Krakow, Poland. 



Abstract 



We compare predictions based on small x (QCD) dynamics with recent data for deep inelastic 
events containing forward jets or forward ir° mesons. We quantify the effect of imposing the (higher 
order) consistency condition on the BFKL equation and study uncertainties inherent in the QCD 
predictions. We also estimate the cross-section for the forward production of two jets. 



1 Introduction 



The positron-proton collider at HERA in the DESY facility has opened a window on a rich vein of 
fascinating physics. In recent years the range of the available kinematic space for deep inelastic pro- 
cesses has been extended ever further, and for a greater breadth of processes, allowing for increasingly 
stringent tests on the physics of Quantum Chromodynamics (QCD), the theory of strong interactions, 
at perturbative energies. 

Pertubative QCD allows us to predict the evolution in kinematic space (but not the initial form) of 
the parton distribution functions. These in turn drive the equations that describe physical observables. 
In distinct regions of (x, Q 2 ) space we have two different modes of parton evolution. At high Q 2 pQCD 
requires we resum the contributions of as log(Q 2 /<5o) t erms - This yields the well documented DGLAP 
equations. As centre of mass energy ^/s, increases at moderate Q 2 , we can technically encounter a 
second large logarithm, this time in 1/x ~ s/Q 2 . Resummation of this type leads to the celebrated 
Balitskij-Fadin-Kuraev-Lipatov (BFKL) equation for the gluon. The BFKL equation corresponds to 
ladder diagrams with (reggeized) gluon exchange along the chain. 

One of the main predictions of the BFKL formalism is singular power law small x behaviour, 
of deep inelastic scattering. The exponent A has been calculated at leading order || and next- 
to-leading order || ^ [j| which contributes large negative corrections to the LO value. The leading 
order approximation of the BFKL equation should therefore be regarded as unreliable. It should 
be emphasised that the subleading log(l/2;) effects are so strong that their next-to- leading (NLO) 
approximation alone is entirely unreliable already for reasonably small, yet relevant values, of the 
strong coupling, that is for as > 0.15. The exact form of the sub-leading contribution, resummed 
to all orders, is unfortunately unknown. It may however be possible to pin down the dominant non- 
leading effects of well-defined physical origin, and to perform their exact resummation (that is going 
beyond the unreliable next-to-leading approximation) ]|, ||. In ref. || we have identified one class of 
such sub-leading corrections which followed from imposing a certain kinematical consistency constraint 
on the available phase space of the emitted gluons along the chain. At the next-to-leading level this 
constraint exhausts about 70% of the corrections to the BFKL exponent A and, more important, the 
effects of this constraint can be resummed to all orders. After including the constraint in the (LO) 
BFKL equation, we obtain a formalism which makes it possible to implement the dominant subleading 
effects and to resum them to all orders. In ref. |7| the BFKL equation supplemented by the consistency 
constraint was used in a quantitative analysis of the structure function F2 within a unified BFKL and 
DGLAP scheme. In this paper we wish to study the effects of these modifications on the properties 
of final state processes in deep inelastic scattering (DIS), which are particularly sensitive on BFKL 
dynamics. These are DIS accompanied by a forward jet Q and the associated process of DIS with 
forward 7r°'s Q. The existing calculations of these processes have been performed within a LO BFKL 
framework, and it is clearly essential to improve the analysis by incorporating the effects of the above 
constraint. 

The DIS + forward jet measurement first proposed by Mueller || can be a very useful tool for 
probing the BFKL dynamics where the diffusion of the transverse momentum along the gluon chain 
plays an important role. By choosing the configuration in which k 2 T ~ Q 2 , where k 2 T denotes the jet 
transverse momentum, we eliminate the phase space region of strongly-ordered transverse momenta of 
the gluons along the chain, i.e. the region k 2 T <C k\ ... <C k 2 -C Q 2 , which is the dominant factor driving 
the increase of F% in the double leading log approximation. Moreover in the forward configuration 

S> x, where x is the Bjorken parameter, we have large subenergy available for jet production that 
justifies the use of BFKL dynamics. 
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We can also study DIS + forward 7r°'s |J. This is a more refined version of the DIS + forward jet 
case, in the sense that we no longer have potential ambiguities derived from hadronisation effects or 
jet finding algorithms. The distinguishing features are discussed more fully in section 5. 

We first motivate higher order constraints which should naturally be implemented in a BFKL 
formalism. We then show how observables for DIS + forward jet/V can be calculated from the 
unintegrated gluon distribution, evolved from the virtual photon end of the BFKL ladder [] , unlike 
the hybrid evolution used in ref. ||. It is more convienient to impose the higher order constraints 
in the present form of the evolution. We then compare our calculations with experimental results for 



DIS + forward jet production from HI [10] and ZEUS [11], and show a comparison of our analysis 
with very recent results for 7r° cross-sections obtained by the HI collaboration at HERA |12| |. Finally 
we estimate the cross-section for production of two jets satisfying the forward criteria. 



2 The unintegrated BFKL gluon from fusion 

In this section we shall describe the formalism needed for the theoretical description of the DIS + 
forward jet measurement based on the BFKL equation [13]. The natural quantity driving physi- 
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Figure 1: Diagram showing the 7*5 fusion, including the soft gluon radiation ladder. 



cal processes in DIS in the small Bjorken-x region of (x,Q 2 ) phase-space is the unintegrated gluon 
distribution, (pi 

x] gi (xlQ 2 ,Q 2 ) = J ^fMx],k 2 T ,Q 2 ), (1) 

where gi(x^, Q 2 ,Q 2 ) are the conventional gluon distributions of a photon of virtuality Q 2 and polari- 
sation i = T,L, probed at a scale Q 2 . The parameter denotes the longitudinal momentum fraction 
of the parent virtual photon carried by the gluon and kx denotes its transverse momentum. The gluon 
four momentum k has the following Sudakov decomposition 

k = -xPp + x]q' + k T , (2) 

1 The present form of analysis in which we iterate from the virtual photon will facilitate comparison with the analysis 
of forward jet production in terms of the partonic structure of the virtual photon. 
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where p is the four momentum of the nucleon, and the light-like vector q' is defined by 



q' = q + xp. (3) 

The variable x is the Bjorken parameter conventionally defined as x = Q 2 /(2p.q) and, as usual, 
Q 2 = — q 2 with q being the four momentum carried by the virtual photon. 

The contribution to <j) arising from the virtual photon can be thought of in terms of virtual photon- 
gluon fusion through a quark-box coupled to an equation describing the subsequent evolution of the 
gluon distribution. The quark box driving terms, 4>t an< ^ ^l^' can ^ e evaluated perturbatively 

<$'(xJ,4,Q 2 ) = 



(xl 4, Q 1 ) = J2 jf ' dP J <fW(l - 13? " ^ 



6 <?>, 
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where the denominator functions are 

D x = k 2 t + 13(1 - I3)Q 2 + m 2 v 

D 2 = (bt - k T ) 2 + 0(1 - [3)Q 2 + m\. 

The quark mass m q is set to zero for the light u, d, and s quarks and taken to be 1.4GeV for the charm 
quark. Taking these driving terms we then evolve the unintegrated gluon distributions, (j)T and (f>L, 
using the BFKL equation from the virtual photon end of the BFKL ladder 

(0), 



<P l (xl,k^Q 2 ) = ^>(xlki,Q 2 )+ 



1 dz f°° dk' 2 <f>i{xVz, k' 2 , Q 2 ) - (kixVz, k 2 T , Q 2 ) <t>i(x1/z, k 2 , Q 2 ) 



' - ' 1-7, "T [ I k T - k 'f I (4k'£ + *4) 2 



r 2 h' 2 I It 2 — h' 2 



(5) 



where i = T or L, as = 3as/7r, and x^ is the fraction of the virtual photon's longitudinal momentum 
carried by the gluon. We will show results for infrared cut-off Kq varied within the range 0.5 — lGeV 2 . 

We require strong ordering in (x^)i along the gluon chain 



(xl) < {xl)x « ... « (xl 



in- 



The present approach differs from our earlier work || in that we no longer have a free parameter zq 
which specifies the end of the BFKL evolution. In our previous treatment |J we performed a hybrid 
evolution of the gluon ladder in the sense that, while we allowed the k\ diffusion from the quark box, 
the longitudinal momentum fractions were defined with respect to the proton end of the ladder and 
strongly ordered along the gluon chain. In that case zq was adjusted to give the correct normalisation 
for the DIS + jet data. In the present calculation the normalisation is, in principle, fixed by the 
theory. However, in practice, there is freedom due to the choice of the value of QCD scale, which we 
will quantify later. In summary, the present description of the behaviour of the quark-box gluon chain 
system can be considered as a calculation of the unintegrated gluon content of the virtual photon. 
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3 Higher order corrections 



Calculations ||, ||, |B| have shown that the NLO corrections to the BFKL equation are large. This casts 
doubt upon the quantitative predictive power of any LO BFKL approach to small x DIS dynamics. 
Clearly, for phenomenological purposes, a method of including the higher order effects should be 
considered. One origin of subleading effects, which contributes to all orders, is the so called consistency 
constraint (CC) |J. This requires that the virtuality of the emitted gluons along the chain should 
arise predominantly from the transverse components of momentum, in order for the small x (small 
x g) approximation to be valid, that is 

4 < f , (6) 

where we have omitted a factor of (1 — z). 

We can motivate the inclusion of the consistency constraint in the analysis by reference to the 
exponent A governing small x behaviour, x~ x . It has been shown that by truncating the all order 
BFKL + CC solution at NLO, one recovers some 70% of the full, explicit NLO calculation in the 
effective exponent. It looks reasonable to suppose that this comprises the dominant part of all higher 
order effects. 

Another physical source of sub-leading contributions comes from the imposition of the conservation 



of energy-momentum in multigluon emission. Such an effect has been investigated in ref. [14|. It turns 
out that the consistency constraint subsumes energy-momentum conservation over a wide region of 
the allowed phase space WA. 
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Figure 2: Diagram to illustrate the radiation of soft gluons from the reggeized gluonic propagator in 
the i-channel. 

The CC places restrictions on the available phase space under the integration in the BFKL equation 
(||). This manifests itself as a ©-function multiplying the kernel component governing real gluon 
emissions. Upon implementing this, we obtain a modified BFKL equation 

Mxl4,Q 2 ) = cpf\x],4,Q 2 )+ 

, , f 1 dz r°° dk' 2 J O (4/k% - z) MxVz, k$, Q 2 ) - <t>i{x1/z, k 2 , Q 2 ) MxVz, k 2 , Q 2 ) 
^ k ''l^l-^< + (4 ,, + 4)5 j 

(7) 



,;. : .//„-' k' 2 ) \k 2 -k' 2 \ (4^ + 4)5 
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4 Forward jet production in DIS. 
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Figure 3: Diagram showing forward jet production driven by 7*3 fusion coupled to evolution through 
a BFKL-type ladder. The struck parton from the proton could in principle be either a gluon or a 
quark. 



Given that measurements of F2 proved unable to unambiguously determine small x dynamics, less 
inclusive alternative approaches were needed. One such process is Mueller's proposal H: 

7* +P — >jet + X 

According to Mueller, DIS events containing an identified forward jet provide a particularly clean 
window on small x dynamics. By requiring that Q 2 ~ k 2 T we ensure minimal DGLAP type evolution, 
and any remaining features must be the result of dynamics from the small x limiting region of the 
(x,Q 2 ) phase-space. Moreover by considering jets with large fractions of the proton's longitudinal 
momentum, x P allowing x/x P to be as small as possible, on the one hand, will allow the use of BFKL 
dynamics, while on the other hand, involve parton distributions at values of x p where they are well 
known. 



4.1 QCD prediction for DIS + forward jet 

The kinematics are such that we may use strong ordering at the parton-gluon-jet vertex 

^aT < ^ ^j'T 1 ' "^j ^3 

where the Sudakov decomposition of the jet four momentum in terms of light-like momenta p and q' 
in analogy to @) is 

kj = x p jP + x]q' + k jT . (8) 

We can relate the longitudinal momentum fractions at the jet vertex through strong-ordering and the 
jet on-shell condition, k 2 = 0. This allows us to deduce that 

2x p jX ](p.q') = k] T , 
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fi2 j,2 ™ 

-v 7 N'T S'T X 

X °- 3 2x p j (p.q) 

Using the prescription described in the previous section, we solve numerically for <j) in the modified 
BFKL equation, via an expansion in Chebyshev polynomials, and finally determine the unintegrated 
gluon distribution 

M^p, tf T , Q 2 ) = U*1 = k] T , Q 2 ). (9) 

x j ^ x j 

&i(x/x P j,k 2 T: Q 2 ) can then be used to calculate the differential structure functions that drive the 

forward jet process through 

(10) 

where J2 a f a is the sum over parton distributions assuming i-channel pole dominance. 

a q 

These quantities can be substituted into an expression for the differential cross-section for forward jets 



dxdQ 2 dx p j dk 2 T xQ 4 
where y = {p.q)/{p e -q) and F 2 = Ft + Fl- 



8 2 F 2 1 2 d 2 F T 



dx v dk 2 T 2 dx^dk 2 T j 



(12) 



4.2 Comparison with HERA data 

Finally we are in a position to be able to confront our calculation with the HERA data. We numerically 
integrate the differential cross-section over the kinematic regions used by the HI [10] and ZEUS 0] 
collaborations, see Table [l]. The results are shown in Fig. 4. Since the parton distributions (|ll|) are 
required in a kinematic domain where they are well known the results are not sensitive to the particular 
set that is used. However there is a dependence of the results on the QCD scale. To be consistent we 
use a recent leading order set of partons. We take the LO set from [O], for which A(QCD) = 0.174GeV 



for four flavours. The three curves in the figure correspond to the arguments of as in (|4j) and (1C) 
and the infrared cut-off in (0) being respectively taken to be ^ 

(i) (k T + m 2 )/4, 4/4, k 2 = 0.5GeV 2 (upper dashed) 

(ii) (k' T + m 2 )/4, k^/A, fcg = lGeV 2 (continuous) (13) 

(iii) (/4 + m 2 ), k T , k 2 = 0.5GeV 2 (lower dashed). 

The sensitivity to the choice of scales for as is therefore seen by comparing the predictions for (i) 
and (iii), and the sensitivity to the value chosen for the infrared cut-off fcg by comparing (i) and (ii). 
We see that the uncertainty due to k 2 is much less than the uncertainty due to the choice of scales. 
From Fig. 4 we see that the shape of the x distributions of the DIS + forward jet data are modelled 
well. Moreover the predicted normalisation is satisfactory in that agreement exists for a physically 
reasonable choice of scales and of the infrared cut-off k^. The experimental data correspond to the 
measurements at the hadron level. The hadronisation effects can lower the cross-section by about 
15-20%. We comment on the comparison with data further in section 7. 

2 When the scale of as is less than in choices (i) and (ii), we freeze the coupling at as(kg). 
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Figure 4: The DIS + forward jet differential cross-section versus Bjorken-a; as measured at the hadron 
level by the HI [^] and ZEUS [11| collaborations. The kinematic cuts are given in Table |l[ The curves 
are predictions at the parton level, based on the BFKL formalism including sub-leading corrections, 
corresponding to the three choices of scales and infrared cut-off given in (\H 
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HI cuts 


ZEUS cuts 


E' e > HGeV 

y e > 0.1 
x P > 0.035 
k jT > 3.5GeV 
0.5 < k] T /Q 2 < 2 
7° < 9 jet < 20° 
160° <6 e < 173° 


E' e > lOGeV 

Ve > 0.1 

x p j > 0.036 
Ejx > 5GeV 
0.5 < E] T /Q 2 < 2 

T]j et < 2.6 



Table 1: Table showing the kinematic restrictions imposed on DIS events at HERA for forward jet 
production. 

5 DIS events containing a forward 7r° 




Figure 5: Schematic showing collinear fragmentation of a forward parton jet of momentum kj into a 
forward pion of momentum zk n . The process is described by fragmentation functions, Di(z,[i 2 ) 

A complementary reaction to Mueller's forward jet process is provided by the production of forward 
7T° mesons in deep inelastic events ||. The driving process here is the same as before, with j*g fusion 
coupled to BFKL-type gluon evolution, ejecting a parton within the proton as a forward jet. However, 
we now allow the parton to evolve into a shower of particles containing a pion collinear to the initial 
parton jet 

7*+p— »ir°+X 

There are advantages in attempting to measure the forward 7r° cross-section as compared to the parent 
forward-jet cross-section. 

• Experimentally one can more unambiguously identify a forward 7r° than a jet. 

• By measuring 7r°'s at relatively low x^^pt-k we effectively collect data for energetic forward jets 
with Xj > x n and kjT > PTn, which might otherwise go undetected. 

• Another consequence of looking for a specific final state is that we eliminate our dependence on 
jet-finding algorithms. We can unambiguously determine the signal due to the measurement of 
a real 7r° instead of worrying about jet definitions. 

• Finally, we eliminate hadronisation uncertainties. All non-perturbative hadronisation effects are 
swept into the fitted parametrisations of the fragmentation functions. 

On the other hand, by requiring a single energetic fragment of the jet we reduce the cross-section. 
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5.1 Including forward jet — > ir° fragmentation 

We use the LO tt^ fragmentation functions of Binnewies et al. They present results of the form 

A(^ 2 ) = A^/xV^a - zf^ 2) (14) 

where i = g,q,q and /i 2 is the fragmentation scale. Here we take the fragmentation scale simply as 
k 2 T . The functions for 7r° are taken to be |(vr + + 7r _ ) distributions 

The form of the fragmentation requires the 7r° to carry z = x w /x p of the jet's longitudinal momen- 
tum, in a direction collinear with the initial parton. With this assumption we can relate the transverse 
momentum of the ir° and the forward jet 

k-nT = zkjT- 

Next we obtain the cross-section for ir° production by convoluting the DIS + forward jet cross-section 
( |l2|) with the 7r° fragmentation parametrisations (14) ||] 



dx^dp\^dxdQ 2 



d(J g „„o, , , 4 



<4 



J Xir J J 

d(J 1 D K ( Z p ) | d °9 nn° ( ,2 N 

dx v 3 dk) T dxdQi q ( ' jT> + drfdkfrdxdQ 2 « [z,KjT) 



\ dx p 3 dk] T dxdQ 2 9 K ' jT) ^ 9 

(15) 

The differential cross-sections on the right-hand-side, daj /dx?dkj T dxdQ 2 with j = q,q,g, are given 
by (|l0|) and (|l^) with the replacement of J2 fa by f q , fq and f g respectively. 

5.2 Comparison with forward ty° data 

There exist two sets of measurement of deep inelastic events containing forward 7r° mesons. First there 
arc the 1994 HI data |0|, and now the recent analysis of the 1996 HI data |l^]. The kinematic cuts 



used to obtain these data samples are summarized in Table |2|. The 7r° spectra are defined through 

1 dn n __ J_5av 
N dx atot dx 

where n n is the number of 7r°'s, and N is the total number of DIS events. The recent HI results [jTJ] 
are fully comprehensive, with histograms showing the differential cross-section for it production as 
functions of x, Q 2 , pt-k and r\^. We compute the ir° cross-section using exactly the same three choices 
of the scales of as and infrared cut-off k 2 , as we used for the DIS + forward jet predictions, see Section 
4.2. 

In Fig. ^ we compare our analysis with the 1994 HI ir° spectra in three x w regimes. We find a 
satisfactory description of the shape, and that the normalisation is best described by choice (iii) of 
the scales of (|l3|). The factor of two overshoot of the data by the previous predictions ||] is gone. The 
reason is the suppresion due to the inclusion of the sub-leading corrections. In Figures 7 through 10 
we compare our calculations of the 7r° differential cross-sections as a function of x, Q 2 , pTn and r] n , 
with the new HI data |12]. As in the DIS + forward jet process, we see that there is good overall 
agreement between the predictions and the DIS + forward 7r° data. 
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1994 HI data 


New HI data 


E' e > 12GeV 
Up > o 1 
x w > 0.01 
PTit > lGeV 
5° < 9 V < 25° 
156° <9 e < 173° 


6 > u„ > 1 

Xn > 0.01 

Pt-k > 2.5GeV 
5° < 9 n < 25° 

2 < Q 2 < 70GeV 2 



Table 2: Table showing the kinematic restrictions imposed on DIS events by the HI collaboration 
for forward 7r° production. For the 1996 data, implicit bounds are imposed on 9 e and E' e by the y e 
restriction. 

The continuous curves, which are the set which best describe the forward jet data, are on average 
some 20% above the forward 7r° data. However the overall agreement is well within the theoretical 
and experimental uncertainties. (Note that no allowance has been made for the uncertainty in the 
fragmentation functions for the pion). Moreover note that the jet data are shown at the hadron level. 
If we were to correct for hadronisation effects, then the parton level jet data are expected to be 15-20% 
lower |17|| , and the overall consistency of the QCD description is even better than shown. 

If we look at the description of the forward tt° data in more detail then we see, from Figs. 7 and 9, 
that the Q 2 behaviour of the DIS + forward tt° data is not reproduced in detail by any single one of 
the three sets of predictions. For instance for Q 2 < 4.5GeV 2 the continuous curve is below the data, 
although in agreement within the errors, whereas for Q 2 > 4.5GeV 2 the curve is above the data by 
about two standard deviations on average. Note that the relative Q 2 behaviour evident between Figs. 
7(b,c,d) simply reflects the comparison of the curves with the data in Fig. 9(a). However, it is worth 
repeating that the overall agreement between the predictions of small x dynamics and the DIS + ir° 
data is much improved since the previous comparison |j. We comment further on the predictions in 
Section 7. 



6 Forward dijet production 



Forward dijet production provides a further complementary measure of small x dynamics ||T 

7* + p — > jeti + jet 2 + X 

The additional second jet is then required to fulfil the same forward experimental cuts as for the simple 
jet data. The process is shown in Fig. 11. We may neglect the effects of soft gluon radiation in the 
rapidity gap, as the experimental constraints are too restrictive to allow significant development of a 



second chain of gluon emission between the jets [18] 

Xj2 ^ Xji ~ 0(1). 

Analogously to the single forward jet case, the dijet cross-section is determined by 



da 



iira 2 



dxdQ 2 dxjidk 2 T1 dxj2dk 2 T2 xQ 4 



(1-2/) 



dF 2 
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+ «r 



dF T 



dxjidk 2 T1 dxj2dk 2 T2 2 dx j\dk 2 Tl dx j2dkj T2 
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Figure 6: The vr° spectra © versus Bjorken-x obtained from 1994 HI data (k|. The curves are 
predictions, based on the BFKL formalism including sub-leading corrections, corresponding to the 
three choices of scales and infrared cut-off given in fll3|). [The restriction x/x w < 0.1 limits the 
comparison to the domain x < 10~ 3 in the upper plot.] 
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Figure 7: The ir° differential cross-section versus Bjorken-x obtained from 1996 HI (preliminary) data 
[JO]. The curves are predictions, based on the BFKL formalism including sub- leading corrections, 
corresponding to the three choices of scales and infrared cut-off given in (]l3| ) . 
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Figure 8: 7r° differential cross-section in transverse momentum, obtained from 1996 HI (preliminary) 
data [12|. The curves are predictions, based on the BFKL formalism including sub-leading corrections, 
corresponding to the three choices of scales and infrared cut-off given in (]l3| ) . 
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Figure 9: The 7r° differential cross-section versus Q 2 obtained from 1996 HI (preliminary) data [p^ |. 
The curves are predictions, based on the BFKL formalism including sub-leading corrections, corre- 
sponding to the three choices of scales and infrared cut-off given in (^) . 
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Figure 10: The tt° differential cross-section versus rapidity obtained from 1996 HI (preliminary) data 
[ p^ l . The upper and lower histograms are the predictions of the BFKL formalism including subleading 
corrections respectively corresponding to choices (ii) and (iii) of the scales and given in ( |l3|) . 
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Figure 11: Diagram illustrating the emission of an extra forward jet in addition to the standard 
Mueller process. 



(17) 



with 



X ^ dx 3l dk 2 T1 dx j2 dk 2 T2 " dW*/** & fcgg^ E /a^-x, *jn). (18) 

Here, <^ is the azimuthal angular separation between the two forward jets, and k 2 is given by 

kl = {k jT i + kjT2) 2 - 

The azimuthal (</>) integration is performed numerically. Experimental jet resolution effects are in- 
cluded by imposing a cut, i? m i n in rapidity-azimuthal angle space: 



(AT,) 2 + (M? > Rmin = I- 

The HI results published in 1994 quote a result of 6.0 ± 0.8 (stat) ±3.2(syst) pb for the total cross- 
section of DIS + 2 forward jet events JHJ. The calculation, including the consistency constraint, 
and imposing all the experimental cuts imposed by the HI collaboration, results in a cross-section of 
5.2,4.8 or 2.7pb depending on the different scales in as in (18) chosen as in ( |l~3| ) respectively. The 



predictions for the two-jet /one-jet ratio give 1.0,1.1 and 0.8% respectively, to be compared with the 
HI measurement of 1.1 ± 0.6%. That is, small x QCD is able to satisfactorily describe the observed 
rate of forward dijet production. 



7 Discussion 

The process of deep inelastic scattering (DIS) with a forward jet has quite a chequered history. The 
original "gold-plated" proposal of using DIS (x, Q 2 ) events containing a forward (xj, k"j T ) jet to study 
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Figure 12: The differential cross-section for (a) DIS + forward jet, (b) DIS + forward 7r° (1996 HI 
preliminary), events. The continuous and dashed curves correspond to the inclusion and omission of 
the consistency constraint. The lower curve of each pair corresponds to choosing the scales of as in 
(||) and ([!(]) to be k\ + mi and respectively. The upper curves have scales 1/4 of these values. In 
all cases the infrared cut-off in (Q) is taken to be kg = 0.5GeV 2 . The jet data are at the hadron level; 



converting to the parton level is expected to lower the jet data points by about 15-20% [17 
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small x dynamics originally dates back to Mueller ||. The idea was to study small x/xj dynamics 
with k 2 T ~ Q 2 so that DGLAP effects are absent. The observation of a forward jet allows the deep 
inelastic scattering to take place off a parton (in a domain where the distribution f a (xj,kj>j) is well 
known) so that the small x/xj dynamics can be studied free from the uncertainties in proton struc- 
ture. The detailed formalism was provided soon after [13]. Then it was the turn of the experiments to 
collect a DIS + forward jet data sample. To obtain sufficient statistics, jets of relatively low k 2 T had 
to be considered, bringing problems of jet identification and measurement. On the theoretical side all 



phenmenological analyses |13|, 19, S] were based on the LO BFKL formalism, with a tendency of the 
predictions to overshoot the observed DIS + forward jet cross-section. At the same time it was also 
shown that one cannot obtain a sufficient increase of the cross-section, with decreasing x, from fixed 
order QCD calculations. In fact the fixed order (NLO) predictions are found to be a factor of 4 or 
more below the data M, M. 



Recently a different approach [£l], 22] to forward jet production has been proposed, based on de- 



scribing the data in terms of the partonic structure of the virtual photon and DGLAP evolution. In 
this picture the dominant contribution in the region Q 2 < kj> „• to the production of forward jets comes 
from the hard scattering of partons in the virtual photon with those in a proton. The calculations 
presented in ref. (2^[ are based on the parton distribution functions of the virtual photon evolved 
from some (low) scale Qq to the scale /i 2 = —t (which is greater than but less than k\j + Q 2 )- 
The parton distributions in virtual photons which are used in the estimates of the forward jet cross- 
sections were taken from ref. p3| . These parton distributions contain a rather arbitrary, although 
phenomenologically plausible, parameterization of the Q 2 dependence. At large Q 2 , however, the 
parton distributions in the virtual photon are dominated by the point-like contribution which can be 
completely specified in perturbative QCD. In fact the approach developed in [21] is not very sensitive 
to the specific parameterisation of the parton distributions of the virtual photon. Our calculations 
based upon the BFKL equation may also be interpreted in terms of hard scattering of partons in a 
virtual photon. In this case, however, the Q 2 dependence of the parton distributions is not arbitrary, 
but rather is dynamically specified by the impact factor defined by the quark box. This impact factor 
corresponds to the point-like j*qq coupling. The analogy with the resolved photon picture becomes 
more evident if we approximate the BFKL evolution by its double leading logarithm approximation 
corresponding to strongly-ordered transverse momenta from Q 2 towards the hard scale kj,j- The 
BFKL-based calculation should not therefore be interpreted as an alternative explanation to the re- 
solved virtual photon picture since the latter is just part of the former. The BFKL description also has 
the merit of treating in a unified way all possible kinematical configurations. It does not in particular 
divide the underlying mechanism depending upon the relations between the potentially large scales 
(Q 2 and fcy ■) of the problem. 



There have been two recent developments, both of which we have considered in this paper. First, 
DIS + forward 7r° process has also been measured. These data provide a complementary measurement 
which overcomes the experimental ambiguities inherent in the measurements of forward jets, albeit 
at the expense of a reduced rate. Second, as far as BFKL theory itself is concerned, there now exist 
complete NLO results available. These indicate that sub-leading effects are very important and cannot 
be neglected in a confrontation of small x dynamics with the data. In this paper we include subleading 
contributions in the description of DIS + forward jet and forward tt° data. The effect of including 
the subleading terms can be seen from Fig. 12. We show the predictions for DIS + forward jet and 
DIS + forward tt° with and without the consistency constraint included for two physically reasonable 
choices of scale. We see that the subleading terms (that is the consistency constraint) reduce the 
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predictions at the smaller values of x by almost a factor of two. There is a sizeable ambiguity in the 
predictions due to the choice of scale, but nevertheless both sets of data favour the inclusion of the 
sub-leading terms. Even though the subleading effect suppresses the cross-sections, the predictions 
remain sufficiently steep with decreasing x to describe the data. By inspecting the comparison shown 
in the two plots in Fig. 12 we note that, relative to the data, the DIS + forward 7T° predictions are 
lower than those for the DIS + jet process. However the comparison for the DIS + tt° process is Q 2 
dependent, see Fig. 9 and 7. On average the forward 7r° predictions are about 20% above the data if 
the scales are chosen so as to give an optimum description of the forward jet data. This discrepancy is 
well within the total uncertainties. Moreover, as was mentioned before, the experimental data on the 
forward jet cross-sections correspond to the measurements at the hadron level, while our theoretical 
predictions concern the parton level. The hadronisation effects are expected to raise the parton level 
jet cross-section by about 15-20% fj"7j| , and so the overall consistency between the forward jet and 7r° 
data and the small x QCD predictions shown at the parton level in Fig. 12, is in fact better than 
shown. Since we can identify and measure less energetic forward 7r°'s than jets, we are able to sample 
smaller values of x in the former process. 

We conclude that there exists an economical, physically-based description of both the DIS + 
forward jet and the DIS + forward 7r° data in terms of the small x QCD framework. 
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